Antibunching photons in a cavity coupled to an optomechanical system 
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We study the photon statistics of a cavity linearly coupled to an optomechanical system via second order 
correlation functions. Our calculations show that the cavity can exhibit strong photon antibunching even when 
optomechanical interaction in the optomechanical system is weak. The cooperation between the weak optome- 
chanical interaction and the destructive interference between different paths for two-photon excitation leads to 
the efficient antibunching effect. Compared with the standard optomechanical system, the coupling between a 
cavity and an optomechanical system provides a method to relax the condition for obtaining single photon by 
optomechanical interaction. 
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I. INTRODUCTION 

Single-photon sources are of considerable interest for im- 
portant applications in quantum information processing (see, 
e.g., review [1] and references therein), for example, key 
distribution in quantum cryptography [2]. Initially, single- 
photons are achieved by attenuated laser pulses with an aver- 
age number of photons per pulse smaller than one. While the 
emitting photons show bunching statistics, and this may result 
in low transmission rates or information leakage in quantum 
key distribution. Clearly, single-photon emitter should deliv- 
ers one, and only one photon at a time, in other words, the 
emitting photons should show antibunching effect. Various 
mechanisms have been employed to arouse antibunching ef- 
fect, and one of them is to utilize photon blockade, which is a 
kind of quantum optical effect of preventing subsequent pho- 
tons from resonantly entering the cavity yfl. Many theoretical 
proposals have been made for generating photon blockade, 
such as, optical cavity with four-level atomic ensemble [3], 
quantum well exciton in a planar microcavity ||4j], coupled 
nonlinear cavities lyHZl], and optomechanical systems [8, 9]. 
To date, photon blockade has been observed in both optical 
cavities flolfTlll and superconducting circuits HHHD]. 

Optomechanical system is attracting more and more atten- 
tion because of many grea t exp erimental achievements in re- 
cent years (see, reviews lll4i [l5ll ). Sideband cooling of a 
mechanical oscillator into its ground state paves the way for 
putting mechanical oscillators into quantum mechanics lfTr3l - 
13011 . The optical response of optomechanical systems to sig- 
nal field can be modified by driven the cavity with a strong 
laser field, leading to effects such as electromagnetically in- 
duced transparency (EIT) ll3ll - l35ll and normal-mode split- 
ting 1 31, 36]. Normal-mode splitting is the evidence for enter- 
ing the strong optomechanical coupling regime. 

In Refs. Jit 0], the statistical properties and blockade of 
photons in optomechanical systems have been theoretically 
studied. These studies show that the nonlinear quantum op- 
tics at single-photon level using optomechanical systems be- 
comes possible under strong optomechanical coupling condi- 
tion. Thus the optomechanical system can act as a single- 
photon source. Moreover, we showed that photon induced 



tunneling can occur in optomechanical systems [37]. Re- 
cently, it was shown that the nonlinear interactions in two cou- 
pled optomechanical systems can be significantly enhanced 
for mechanical frequencies nearly resonant to the optical level 
splitting l38l l39ll . However, the photon blockade effects still 
only appear in the strong coupling regime, which is beyond 
the reach of most experiments in the single photon regime. 
Thus there is a question whether single-photon states can be 
generated using weak optomechanical interaction. 

Recently, Liew and Savona analyzed two coupled nonlin- 
ear cavities |6|. they found that the single-photon statistics 
can exhibit strong antibunching in such coupled-mode sys- 
tems with weak Kerr nonlinearity compared to the single- 
mode case. Later on, such strong antibunching is attributed 
to the destructive quantum interference effect [7], and authors 
further extended their theory to two coupled-cavities with 
a two-level quantum emitter, which is embedded in one of 
cavities. They theoretically demonstrated that perfect photon 
antibunching can be obtained even for single-atom cooperativ- 
ity on the order of or smaller than unity. These studies |]6l 0] 
have opened up a door towards nonlinear quantum optics at 
single-photon level using weak nonlinear coupling. 

Motivated by studies in Refs. J(jl Q] and also recent theo- 
retical and experimental progress in coupled-cavity array and 
optomechanical systems, we now study the properties of the 
light field inside a cavity, which is coupled to an optomechan- 
ical system. The paper is organized as follows. In Sec. II, 
the model Hamiltonian is introduced. In Sec. EH, quantum 
Langevin equations of cavity fields and mechanical mode are 
given, and the statistical properties of the cavity fields are 
studied using the second-order correlation function under the 
semiclassical approximation in Sec. IV. In Sec. V, we further 
analyze the second-order correlation function by numerical 
simulation via the master equation, and compare those with 
the results of the semiclassical approximation. Finally, sum- 
mary and conclusions are given in Sec. VI. 



II. MODEL 

As schematically shown in Fig. Q] the system consists of 
two coupled cavities (A and B) with the coupling constant 
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FIG. 1: Schematic diagram for an optical cavity (cavity A) coupled to 
an optomechanical system which consists of an optical cavity (cavity 
B) with an oscillating mirror at one end. Cavity A is driven by a 
weak coherent laser field. 



go. The cavity can be a transmission line resonator, a toroidal 
microresonator, a cavity with two mirrors, or a defect cavity in 
photonic crystal. Without loss of generality and for simplicity, 
we will focus on the cavity with two mirrors. Cavity A is 
driven by a weak probe field with frequency uj c , and cavity 
B consists of an oscillating mirror at one end, modeled as 
a quantum mechanical harmonic oscillator. In other words, 
we study a coupled system, which consists of a driven cavity 
and an optomechanical system. The Hamiltonian of the whole 
system in the rotating wave approximation is given as 

H = hui a a) a + hujbb^b 



+hg (a^b + tfa) + hGb^b\ 
+ihe c (a)e- lu ^ - ae 1 ^) , 



(1) 



where a (a^) is the annihilation (creation) operator for the light 
mode of the cavity A with frequency uj a , b (fet) is the annihi- 
lation (creation) operator for the light mode of the cavity B 
with frequency cjj,, and c (c^) is phonon annihilation (creation) 
operator of the mechanically vibrational mode for the mirror 
with frequency ui m . The parameter G denotes the coupling 
strength between the cavity B and the oscillating mirror, and 
e c presents the coupling strength between the driven field and 
cavity field inside the cavity A. As u> a ~ u>b ^> u> m ,go, we 
have dropped the rapidly varying terms (ab and a'w) corre- 
sponding to the rotating wave approximation. 

To remove the time-dependent factor, let us transform 
the Hamiltonian in Eq. ([TJ into the rotating reference frame 
through a unitary operator R(t) — exp[— iuj c t{(va + b'b)], 
and thus the Hamiltonian in Eq. (fl} becomes 



l c t c 



H = KA a a)a + HA b tfb + Huj, 

+hg (a^b + tfa) +hGtfb(c 
+ihe c (a/ — a) , 



(2) 



where A a = u a — u c and A 5 = w& — uj c are the detunings of 
the frequencies of cavity fields from that of the driving field. 



III. LANGEVIN EQUATIONS AND SEMICLASSICAL 
APPROXIMATION 



The dynamics of the cavity fields and mechanical oscillator 
can be described by quantum Langevin equations. By con- 



sidering the dissipation and fluctuation of the light fields of 
two cavities and mechanical mode, we can write out a set of 
nonlinear quantum Langevin equations as follows 
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where j a , jf, and j m are the damping rates of cavity A, cavity 
B, and the moving mirror, respectively, q = (c + c^), 

p = r^= (c — c^), and gi, = \J~2G. £ is a Brownian stochastic 

force with zero mean value, i.e. (£ (i)) = , which comes 
from the coupling of the oscillating mechanical resonator to 
its thermal environment and satisfies correlation [40-43] 
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where fc^ is the Boltzmann constant and T is the effective 
temperature of the environment of the mechanical resonator. 
<2i n and 6; n represent the vacuum radiation noises input to the 
cavity A and B with (a- m (t)) = (6; n (t)) = 0, and they obey 
the following correlation functions [44] 
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(9) 
(10) 
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Here, we have assumed that whole system is in a low tempera- 
ture environment, and therefore the equilibrium mean thermal 
photon numbers in two cavities at optical frequencies have 
been neglected. 

The dynamic of the system is usually determined by the 
small fluctuations when the system reaches the steady-state. 
Thus, let us now apply semiclassical approximation to solve 
the steady-state with small quantum fluctuations. That is, we 
assume a = a + da, b = f3 + 5b, q = q + Sq, here an, (3q 
and qo are the mean values of the cavity fields and mechanical 
mode when the system reaches steady-states, and operators 
8a, 6b and Sq describe the small fluctuations around steady- 
state and satisfy the condition (Sa) = 0, (Sb) = and (Sq) = 
0. The steady-state values satisfy the following equations 



(y + «A a ) an + igoPo = £c 
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Here, we have used the factorization assumption, e.g., (qb) 
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(q)(b). The dynamics of small fluctuations around steady - 
state can be described by linearizing Eqs. (|3]|6) as 
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here, the high order terms of small fluctuations, e.g., <5g<5a, 
have been neglected. The system is stable only if all the eigen- 
values of the coefficient matrix of the above differential equa- 
tions have negative real parts, and the stability condition can 
be given explicitly by using the Routh-Hurwitz criterion 114 511 . 
However, it is too cumbersome to be given here. All the pa- 
rameters we will use satisfy the stability condition, and it is 
easy to fulfill it for the driven field in our system is weak. 

By applying the Fourier transform and solving dynamical 
equations in the frequency domain, we can obtain 

Sa(oj) = E (u) a in (u) + F (to) a\ n (w) 
+G (w) b in (w) + H (w) b\ n («) 
+Q («) £ (w) , 
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here, we introduce A 6 = A& + <?bgo — g 2 |/So 
dynamical response function of the mirror 
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IV. SECOND-ORDER CORRELATION FUNCTIONS AND 
STATISTICAL PROPERTIES OF CAVITY FIELD 

In this section, we are going to investigate photon statistics 
of cavity A by using the analytical expression given in the last 
section. By taking a = ao + da, the equal-time second-order 
correlation function g aa ' (0) of the light field in the cavity A 
can be given as 



9 (2) 
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where the third term comes from the four-operator correlation 
by using properties of the Gaussian process ifjol RUl . and the 
second term is only one which may take negative value. Using 
the expression of Sa(t), Eqs. ( fT9l . and taking the correlations 
Eqs. d7lfTTI). the correlation of 5a(t) and 8a' (i) are given as 
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For comparison, before analyzing the statistical properties 
of light field in cavity A, we show the equal-time second- 
order correlation function g( 2 > (0) of the standard optome- 
chanical system in Fig. |2ja). The standard optomechanical 
system consists of a cavity with an oscillating mirror at the 
one end SUt]. Figure [2 a) shows that strong antibunching 
appears in standard optomechanical system under strong op- 
tomechanical coupling condition G > 7, and this agrees with 
Ref.Hl. 
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FIG. 2: (a) Equal-time second order correlation function of standard 
optomechanical system gr*- 2 ' (0) plotted as functions of A and G. (b) 
Equal-time second order correlation function of cavity A, gaS (0) 
given by Eq. <3U , plotted as functions of A and G for go = 37. 
The other parameters are 7 a = jt — 7, £c = O.OI7, ui m — IO7, 
to m /-y m = 1000, and T ~ K. 
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FIG. 4: (a) Equal-time second order correlation function of cavity A, 
Qal (0), plotted as functions of G for A = — O.247. (b) Dependence 
of g { aJ (0) on A for G = O.67. The solid curves are calculated 
by master equation, given by Eq. d38| l, and the dashed curves are 
obtained from quantum Langevin equations, given by Eq. d31t . The 
other parameters are the same as in Fig. |2](b). 




1 G/y 



FIG. 3: Equal-time second order correlation function of cavity A, 
gaa (0) given by Eq. d31t . plotted as functions of both g a and G for 
A = — O.247. The other parameters are the same as in Fig. |2](b). 



Now, let us focus on the statistical properties of light field 
in cavity A. From now on, we assume that the two cav- 
ities (cavity A and B) are the same as each other, that is, 
-f a = Jb = 1, w a = w u , the detunings are denoted by 
A a = A;, = A, and normalize all the parameters to 7. The 
equal-time second-order correlation function of the photons 
inside cavity A, g a 2 J (0), is given in Fig.HJb). We observe that 
strong antibunching appears near the point for A = — 0.247 
and G = O.67 in cavity A as go = 37. That is, the photons in 
the cavity can exhibit strong antibunching when it is coupled 
to an optomechanical system under weak condition G < 7. 

The optomechanical interaction in cavity B and the quan- 
tum interference effect between the two cavities (cavity A 
and cavity B) are responsible for the photon antibunching ef- 
fect 0| , The interference is between two paths for two-photon 
excitation in cavity A: (a) the direct excitation from one pho- 
ton to two photons inside the cavity; and (b) one photon tun- 
neling from cavity A to cavity B, then exciting another photon 
in cavity A, and finally the photon inside cavity B tunneling 
back to cavity A. Thus the destructive interference between 
two paths reduces the probability of two-photon excitation in 
cavity A. 

Two-dimensional plot of the equal-time second-order cor- 

(2) 

relation function of cavity A, g aa (0), as functions of both go 
and G is shown in Fig.[3]for A = — O.247. With the increas- 




FIG. 5: Equal-time second order correlation function of of cavity A, 
gaa (0) , plotted as functions of A for different temperatures at G = 
O.67 and 7/27T = 1 MHz. (a) The results are obtained by master 
equation approach, given by Eq. d38t , (b) the results are obtained by 
quantum Langevin equations method, given by Eq. d31t . The other 
parameters are the same as in Fig.|2](b). 



ing of g , the value of G for getting the strong antibunching 

(2) 

gaa (0) is drawn back from the strong optomechanical interac- 
tion condition to the weak optomechanical coupling domain. 
This implies that destructive quantum interference effect in 
the coupled system can be used to lower the strength of the 
optomechanical interaction that is required to achieve strong 
antibunching. 



V. NUMERICAL SOLUTION BY MASTER EQUATION 

In this part, we calculate the second-order correlation func- 
tion by numerically solving the master equation of the density 
matrix, and compare the results to the predictions given by 
analytical solution derived above. The master equation of the 
coupled system is given as [46] 
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where n m is the mean thermal phonon number of the mov- 
ing mirror given by the Bose-Einstein statistics n m = 



5 



(a) 



- ij=aa 
■ ij=ab 

- ij=bb 



(b) 

ij=ab 

ij=bb 









G/y 



A/y 



(21 

FIG. 6: (a) Equal-time second order correlation function (0) 
plotted as functions of G at A = —0.247, ij = aa for cavity A 
(solid), ij = ab for the cross correlation between the two cavities 
(dashed), and ij = bb for cavity B (dashed-dot). (b) Dependence of 
sfj (0) on A for G = O.67. The other parameters are the same as 
inFig.|2](b). 



[exp(Su m / ksT) — 1] _1 , fcg is the Boltzmann constant and T 
is the effective temperature of the moving mirror. The master 
equation can be solved in the basis of the photon and phonon 
number states \n a ,rib,n m ), and p can be written as density 
matrix 

(t)\n a , rib, n m )(n' a , n' b , n' m \. (37) 

If the elements of the steady state density matrix, 
Pn a ,n b ,n m ;n< a ,n' b ,n' m , are given, the equal-time second-order 
correlation function can be easily calculated as 



ffS (0) 
*2? (0) 



Tr[pa* a ] 
[Tr(pa^a)} 2 ' 
Tr[ptf 2 b 2 } 
[Tr(ptfb)] 2 ' 

Trfjoa^fea] 
[Tr( / oat a )][Tr( i ofct&)]' 



(38) 
(39) 
(40) 



(2) 

where g ab (0) is the cross correlation between the photons in 
cavity A and B. 

For comparison, the second order correlation functions cal- 
culated by the master equation and quantum Langevin equa- 
tions are shown in the same figure as functions of G for 
A = — 0.247 in Fig.|4](a), and as functions of A for G = O.67 
in Fig. [4] (b). From Fig. (a), we can see that the results 
obtained by the two methods agree with each other as G is 
small, especially for G < O.67. As shown in Fig. 0(b), for 
G = O.67, the predictions by the two methods match quanti- 
tatively. But with the further increasing of G, the difference 
between them becomes significant, and linearized quantum 
Langevin equations method can only describe this qualita- 
tively. 

The second order correlation functions calculated by the 
master equation and quantum Langevin equations for differ- 
ent temperatures are shown in Fig. From Fig. 0we can see 
that the increase of the phonons will suppress the exhibition of 



antibunching effect. In order to get strong antibunching, the 
mean phonon number of the environment must be small be- 
cause the phonons in the environment may disturb the quan- 
tum statistics of the system. Compare the results in Fig. 0(a) 
and (b), in low temperature (T = 0.1 mK), the results ob- 
tained by the two methods are close; as the temperature goes 
higher (T = 2 mK), the discrepancy between them becomes 
larger. 

Finally, let us give a glance of the statistical properties of 
photons in the entire system. The equal-time second order cor- 
relation functions gj9 (0) can be calculated by using Eqs. (|38l 
l40l and the results are shown in Fig. From Fig. 0(a) we 
can see that, under weak optomechanical interaction condi- 
tion (G < 7), there is strong antibunching in cavity A around 
G = O.67, while weak antibunching in cavity B and bunching 
for the photons between the two cavities. 

Dependence of g^ (0) on A is drawn in Fig. 0(b). Fig.0 
(b) shows us two interesting phenomena as A is in differ- 
ent domains: For A < — O.O87, there is strong antibunch- 
ing in cavity A and weak antibunching in cavity B, while the 
cross correlation between the modes in the two cavities ex- 
hibits bunching. On the contrary, when A > O.I7, there is 
bunching in cavity A and cavity B, while the cross correla- 
tion between the modes in the two cavities exhibits weak an- 
tibunching .9^(0) < 1. Under weakly driven condition ll37ll . 



9aa(0) < 1 and g^'(0) < 1 indicate that there is no more 

(2) 

than one photon in each cavity, and g ab {0) > 1 shows that 
there is big chance that each cavity has one photon simulta- 

(2) (2) 

neously. g aa (0) > 1 and gL, (0) > 1 indicate that there is 

big chance for more than one photon present in each cavity, 

(2) 

while g ab (0) < 1 shows that the probability that each cavity 
has one photon simultaneously is low. In other words, when 
the system is driven weakly, and there are two photons in the 
coupled system, if A < — O.O87, they are likely to be in the 
state that each cavity has one photon simultaneously, and if 
A > O.I7, they prefer to stay in one of the cavities together at 
the same time. 



,(2) 



VI. CONCLUSIONS 

We have studied the photon statistics of a cavity linearly 
coupled to an optomechanical system. Due to destructive 
quantum interference effect between different paths for two- 
photon excitation, the cavity can exhibit strong photon an- 
tibunching with weak optomechanical interaction in the op- 
tomechanical system. Both analytical and numerical methods 
are employed to figure out our result. The result brings hope 
to us of observing blockade effect with current experimental 
parameters of optomechanics. 

We gratefully thank Yu-xi Liu for comments. 
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